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Abstract 

Motivated by the recent indications of the possibility of sizable deviations of the mixing- 
induced CP violation parameter, Sf, in the penguin-dominated b — ► sqq transition decays such 
as B° — > ((/), lo, p , if, 7], 7T°, fo)Ks from sin 2/3 determined from B — ► J/ipK$, we study final-state 
rescattering effects on their decay rates and CP violation. Recent observations of large direct 
CP asymmetry in modes such as B° — » K + ir~ , p~ir + means that final state phases in 2-body 
B decays may not be small. It is therefore important to examine these long-distance effects on 
penguin-dominated decays. Such long-distance effects on Sf are found to be generally small [i.e. 
0(1 — 2%)] or negligible except for the loKs and p°K$ modes where Sf is lowered by around 15% 
for the former and increased by the same percentage for the latter. However, final-state rescattering 
can enhance the ojKg, 4>Ks, rj'Ks, p°K$ and ir°Ks rates significantly and flip the signs of direct 
CP asymmetries of the last two modes. Direct CP asymmetries in uKg and p°K$ channels are 
predicted to be AwK s ~ —0.13 and A p o Ks « 0.47, respectively. However, direct CP asymmetry 
in all the other b —> s penguin dominated modes that we study is found to be rather small ( < a 
few percents), rendering these modes a viable place to search for the CP-odd phases beyond the 
standard model. Since ASf(= —rjfSf — Sj/^ Ks , with r/f being the CP eigenvalue of the final state 
/) and Af are closely related, the theoretical uncertainties in the mixing induced parameter Sf and 
the direct CP asymmetry parameter Af are also coupled. Based on this work, it seems difficult to 
accommodate |A>S/| > 0.10 within the SM for B° — > (0, u, p°, ?/, rj, tt°)Ks, in particular, rj'Ks is 
especially clean in our picture. For foK$, at present we cannot make reliable estimates. The sign 
of the central value of ASy for all the modes we study is positive but quite small, compared to the 
theoretical uncertainties, so that at present conclusive statements on the sign are difficult to make. 



I. INTRODUCTION AND MOTIVATION 



Possible New Physics beyond the Standard Model is being intensively searched via the mea- 
surements of time-dependent CP asymmetries in neutral B meson decays into final CP eigenstates 
defined by 

; — ; ; : ; — r = 6 f Sill Ami) + Af COS! Ami), (1.1) 

where Am is the mass difference of the two neutral B eigenstates, Sf monitors mixing-induced 
CP asymmetry and Af measures direct CP violation (in terms of the BaBar notation, Cf = —Af). 
The CP- violating parameters Af and Sf can be expressed as 



1-|A/| 2 _ 21mA/ 

Af ~~i + \x^' Sf -TTJXfT^ (L - ) 

where 



_ q B A(B° -> f) 
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In the standard model A/ r\fe~ 2% ^ [see Eq. (j2.14j) below] for 6 — ► s penguin-dominated or pure 
penguin modes with rjf = 1 (—1) for final CP-even (odd) states. Therefore, it is expected in the 
Standard Model that —rjfSf ~ sin 2(3 and ~ with (3 being one of the angles of the unitarity 
triangle. 

The mixing-induced CP violation in B decays has been already observed in the golden mode 
P° _ mKa for seV eral years. The current average of BaBar .1] and Belle 3 measurements is 

sin 2(3 « 5"j /v ,^ s = 0.726 ± 0.037 . (1.4) 

However, the time-dependent CP-asymmetries in the b — ► sgg induced two-body decays such as 
P° — ► (0, 7T°, r/, fo)Ks are found to show some indications of deviations from the expectation of 
the Standi Mode, (SM). The BaBar B and Belle Q resuhs and then averages are shown in Table 
J] In the SM, CP asymmetry in all above-mentioned modes should be equal to Sj/^x with a small 
deviation utmost O (0.1) 5]. A S discussed in g. this may originate from the 0(A 2 ) truncation 
and from the subdominant (color-suppressed) tree contribution to these processes. From Tabled 
we see some possibly sizable deviations from the SM, especially in the ij'Ks mode in which the 
discrepancy AS^Ks = —0.30 ± 0.11 is a 2.7a effect where 

ASf = -r]fSf - Sj/^ Ks . (1.5) 

If this deviation from Sj/^x is confirmed and established in the future, it may imply some New 
Physics beyond the SM. 

In order to detect the signal of New Physics unambiguously in the penguin b —* s modes, it is 
of great importance to examine how much of the deviation of Sf from Sj/^k is allowed in the SM 



ot great importai 

BEBBBQ 



]. In all these previous studies and estimates of ASf, effects of FSI were not taken 
into account. In view of the striking observation of large direct CP violation in B° — ► K^ir^, it is 
clear that final-state phases in two-body B decays may not be small. It is therefore important to 
understand their effects on ASf. 
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TABLE I: Mixing-induced and direct CP asymmetries —rjfSf (first entry) and Af (second entry), 
respectively, for various penguin-dominated modes with rjf being the CP eigenvalue of the final 
state. Experimental results are taken from ^, 



Trivial St.at.p 




BaBar 131 


Belle 141 


J. 1 V KsL. L i t-, v> 






50 ± 25 + S'2T 


08 ± 33 ± 09 


35 ± 20 




— U 


00 ± 0.23 ± 0.05 


0.08 ± 0.22 ± 0.09 


0.04 ± 0.17 


uK s 




0.50±g;|| ± 0.02 


0.76 ± 0.65igi| 


55+ - 30 
u - oo -0.32 






0.56±g;|? ± 0.03 


0.27 ±0.48 ±0.15 


0.48 ±0.25 







30 ±0.14 ±0.02 


0.65 ±0.18 ±0.04 


0.43 ±0.11 







21 ±0.10 ±0.02 


-0.19 ±0.11 ±0.05 


0.04 ± 0.08 


ir°K s 




0.35±g;|§ ± 0.04 


0.32 ±0.61 ±0.13 


34+ - 27 




-0 


06 ±0.18 ±0.03 


-0.11 ±0.20 ±0.09 


-0.08 ±0.14 


foK s 




0.95±g;|l ±0.10 


-0.47 ±0.41 ±0.08 


0.39 ±0.26 







24 ±0.31 ±0.15 


-0.39 ±0.27 ±0.09 


-0.14 ±0.22 



The decay amplitude for the pure penguin or penguin-dominated charmless B decay in general 
has the form 

M(B° - /) = V ub V: s F u + V cb V c * s F c + VthV^F 1 . (1.6) 

Unitarity of the CKM matrix elements leads to 

M(B° - /) = V ub V: s A u f + V cb V c * s A c f « AX 4 R b e- l ^A u f + AX 2 A}, (1.7) 

where A u f = F u - F t , A% = F c - F*, R b = \V ud V ub / {V cd V cb )\ = J? + W with p, fj, A, A being the 
Wolfenstein parameters [ll|. The first term is suppressed by a factor of A 2 relative to the second 
term. For a pure penguin decay such as B° — > <pK°, it is naively expected that AJ is in general 
comparable to Aj in magnitude. Therefore, to a good approximation —rjfSf ~ sin 2/3 Sj/^k. 
For penguin-dominated modes such as ojK$, p°Ks,tt°Ks, Aj also receives tree contributions from 
the b — > uus tree operators. Since the Wilson coefficient for the penguin operator is smaller than 
the one for the tree operator, Aj could be significantly larger than A C f. As the first term carries 
a weak phase 7, it is possible that Sf is subject to a significant "tree pollution". To quantify the 
deviation, it is known that to the first order in rj = (\ u AJ)/(X c Af) pi. Il^ 

ASf = 2\rf\ cos 2/3 sin 7 cos 5/ , Af = 2|rj| sin 7 sin Sf, (1-8) 

with 5f = arg(Aj/A c j). Hence, the magnitude of the CP asymmetry difference ASf and direct 
CP violation are both governed by the size of Aj/A c j. However, for the aforementioned penguin- 
dominated modes, the tree contribution is color suppressed and hence in practice the deviation of 
8, is expected to be small 0. (However , we shall see below that a sizable ASf can occur in ujKs 
and p°Ks modes. For a review of model calculations of ASf, see 010 

Since the penguin loop contributions are sensitive to high virtuality, New Physics beyond the 
SM may contribute to Sf through the heavy particles in the loops (for a review of the New Physics 
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sources contributing to Sf, see ^J). Another possibility is that final-state interactions are the 
possible tree pollution sources to Sf. Both AJ and Aj will receive long-distance tree and penguin 
contributions from rescattering of some intermediate states. In particular there may be some 
dynamical enhancement of light u-quark loop. If tree contributions to AJ are sizable, then final- 
state rescattering will have the potential of pushing Sf away from the naive expectation. Take the 
penguin-dominated decay 5° — > ujK° as an illustration. It can proceed through the weak decay 
B° — > K*~ir + followed by the rescattering K*~tt + — ► ujK° . The tree contribution to B° — ► ET*~7r + , 
which is color allowed, turns out to be comparable to the penguin one because of the absence of the 
chiral enhancement characterized by the a§ penguin term. Consequently, even within the framework 
of the SM, final-state rescattering may provide a mechanism of tree pollution to Sf. By the same 
token, we note that although B° — ► (pK° is a pure penguin process at short distances, it does 
receive tree contributions via long-distance rescattering. 

In this work, we shall study the effects of final-state interactions (FSI) on the time-dependent 
CP asymmetries Sf and Af. In we have studied the final-state rescattering effects on the 
hadronic B decays and examined their impact on direct CP violation. The direct CP-violating 
partial rate asymmetries in charmless B decays to tttt/ttK and pir are significantly affected by final- 
state rescattering and their signs are gene rally different from that predicted by the short-distance 



approach such as QCD factorization [lfiL I18L Il9l |. Evidence of direct CP violation in the decay 
B° — > K~ir + is now established, while the combined BaBar and Belle measurements of B° — ► p^ir^ 
imply a 3.6a direct CP asymmetry in the p + ir~ mode [l7]. In fact, direct CP asymmetries in these 
channels are a lot bigger than expectations (of many people) and may be indicative of appreciable 
LD rescattering effects, in general, in B decays. Our predictions for CP violation agree with 
experiment in both magnitude and sign, whereas the QCD factorization predictions (especially for 
P + tt~) seem to have some difficulties with the data. 

Besides some significant final-state rescattering effects on direct CP violation, another motiva- 
tion for including FSIs is that there are consistently 2 to 3 a deviations between the central values 
of the QCDF predictions for penguin-dominated modes such as B — > K*tt, K<f>, K*<p, Krj' and the 
experimental data This discrepancy between theory and experiment for branching ratios may 
indicate the importance of subleading power corrections such as FSI effects and/or the annihilation 
topology. 

Since direct CP violation in charmless B decays can be significantly affected by final-state 
rescattering, it is clearly important to try to take into account the FSI effect on the mixing-induced 
and direct CP asymmetries Sf and Af of these penguin dominated modes. The layout of the 
present paper is as follows. In Sec. II we discuss the short-distance contributions to the b — > sqq 
transition decays B° — ► (0, u, p°, 7r°, 7/, rj, fa)K$ within the framework of QCD factorization. We 
then proceed to study the final-state rescattering effects on CP asymmetries Sf and Af in Sec. 
III. Sec. IV contains our conclusions. 
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II. SHORT-DISTANCE CONTRIBUTIONS 



A. Decay amplitudes in QCD factorization 



We shall use the QCD factorization approach [16|, [18|, [19( to study the short-distance con- 
tributions to the decays i?° —* ((f), u, p°, 7r°, rj', rj, fo)K . In QCD factorization, the factorization 
amplitudes of above-mentioned decays are given by 



(<PK°\H cS \B°) = ^{ PB -e%) £ A p < 

V ^ p=u,c 



(a 3 + a B )(K°<t>) + « + 4al)(K\) -\(a 7 + a 9 )(R°^ 



x 2f <t) m <t> F* K (ml) + f B f K U 



bs(K°<f>) - lb™(K°c 



2m„ 



m 



2 ( ' 



p=u,c 



a 2 {K Q u)5l + 2(a 3 + a 5 )( K°u) + \{a 7 + a 9 )(K°oo)] 



(p°K°\H cS \B ) 



x 2f w m w F? K (mi) + 

x 2f K m UJ A^{m 2 K ) + f B f K f l 
Gf 



2m u 1 



& 3 (^V^f W (^°) 



2 f ' 



p=u,c 



a 2 (K%)5Z + ^(a 7 + a 9 )(K°p)] 



x 2f p m p F* K (m 2 p ) 



x 2f K m p AQ p (m 2 K ) - f B fKfp 



2m n 1 



m 



2 ( ' 
B J 



(7r ^°|ff eff |5 u 



p=«,c 



a 2 (^V)^ + ^a 3 ,Ew(^V) 



.Gf 



-a 



4,EW 



(7T Z°) 



/^(m|,)(m|-m2) 



a 2 (^V)^ + 2a 3 (KV; 



+ ^3,Ew(if ^ 



a 2 (^)<5? + a 3 (^V) 



+ # 



a3(K°r]' s ) + al(K°r]' s ) - ^a^ E w (K° i]' s ) - \oP A;EIN {rj[ 8 K°) J 



+ (m 2 K )(m 2 B - ml,) 



f_K_ 

V2 



a p A {v'K°) - ra^^) 



+ f B f K (—f«+% 



V2* 



b^'K )-^^ ) 



5 



{f K°\H cS \B c 



Gf 

7i 



p=u,c 



r«a p 6 )(f K°) 



10 



r?o8)(/o* 



x fKF^im^ml 



fsfxifs + fu) 



+ (2a; 
1 



a p 8 )(K"f )f s 



Fo BK (m%)(ml 



WoK") - ^-bf w (f K 



(2.1) 



where Fi t Q and Aq denote pseudoscalar and vector form factors in the standard convention |2 

A p = VpbV* s , 

a 3 {M 1 M 2 ) 
«3,ew(MiM 2 ) = 
and 



a 3 (M 1 M 2 )-a 5 (M 1 M 2 ), 
-o 7 (MiM 2 ) +a 9 (MiM 2 ), 



o? A {M x M 2 ) 



<(MiM 2 ) + rf 2 o^(MiM 2 ), 



(2.2) 



a^ EW (M!M 2 ) = a? (MiM 2 ) + r^a^M,) 



rtip) 



2mi 



2m^ 



4 {H) 



2m 2 , 



1 



f 9 



(2.3) 



m b (fj,)[m s (n) + m q (fj,)]' x 2m h {n)m s {n) 
with the scale dependent transverse decay constant fy being defined as 

(V(p,e*)\q^vq'\0) = f^ip^l - Pl/ e*). (2.4) 

— — 

Note that the a6 penguin term appears in the decay amplitude of B — > 4>K owing to the non- 
vanishing transverse decay constant of the <f> meson. The scalar decay constant f q of /o(980) in 
Eq. (|2.1I) is defined by {fo\qq\0) = mf f q . The decay amplitude for 7J° — ► 7]K° is obtained from 
B° — » 7]'K° by replacing rf — » 77 . Note that the use of nonzero g 2 in the argument of form factors 
in Eq. (|2.1|) means that some corrections quadratic in the light quark masses are automatically 
incorporated. 

The effective parameters a% with p = u,c can be calculated in the QCD factorization approach 



|16j |. They are basically the Wilson coefficients in conjunction with short-distance nonfactorizable 
corrections such as vertex corrections and hard spectator interactions. In general, they have the 
expressions 0, 0] 



af(MiM 2 ) 



4vr 2 

Vi(M 2 ) + _^(MiM 2 ) 



(2.5) 



where i = 1, • • • , 10, the upper (lower) signs apply when i is odd (even), q are the Wilson coefficients, 
C F = (JV 2 - 1)/(2JV C ) with N c = 3, M 2 is the emitted meson and M\ shares the same spectator 
quark with the B meson. The quantities Vi{M 2 ) account for vertex corrections, Hi(M\M 2 ) for hard 
spectator interactions with a hard gluon exchange between the emitted meson and the spectator 
quark of the B meson and Pi(M 2 ) for penguin contractions. The explicit expressions of these 
quantities together with the annihilation quantities 63 and 63 w can be found in }lfil. Il^|. 

For B — » 77' - 'K decay, rj q and rjl in Eq. (|2.1|) refer to the non-strange and strange quark states, 
respectively, of rf'. The decay constants f q '^ c are defined by 

■j=f* ( >)Pii, (V { '\p)\s7^5s\0) = -iffap?, (v { 'Hp) I c7m75c|0) = -if^Pp, 

(2.6) 



fa (,) (p) 1 97^7591 0) 
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with q = u or d. Numerically, we shall use 21 ] 

f q q , = 89 MeV, f°, = 131 MeV, 
/« = 108MeV, /* = -HlMeV, 



./), 



-6.3 MeV, 
-2.4 MeV, 



(2.7) 



recalling that, in our convention, f n = 132 MeV. As for the B — > rf^ form factor, we shall follow 



18] to use the relation F, 







(f 9 (')/ f S ('))^ n t° obtain the values of the form factors F 



BrjO 



The wave functions of the physical r/ and n states are related to that of the SU(3) singlet state % 
and octet state rjs by 



?78 sin <j> + t]q cos i 



n = r/g cos « 



r]o sin < 



(2.8 



with the mixing angle (f> = —(15.4 ± 1.0)° |21. j 

The decay B — > /o(980)if has been discussed in detail in ;22|. Since the scalar meson /q(980) 
cannot be produced via the vector current owing to charge conjugation invariance or conservation of 
vector current, the tree contribution to B — ► foK vanishes under the factorization approximation. 
Just as the SU(3)-singlet r/o, /o(980) in the 2-quark picture also contains strange and non-strange 
quark content 



l/o(980)> 



\SS) COS I 



inn) sin ( 



(2.9) 



with nn = (uu + dd) / y2. Experimental implications for the mixing angle 6 have been discussed in 



detail in 



22 



I; it lies in the ranges of 25° < 6 < 40° and 140° < 6 < 165°. Based on the QCD 
sum-rule technique, the decay constants f s and f n defined by (/o|^|0) = mf f q with /q 1 = nn and 
/o = ss have been estimated in [i^l by taking into account their scale dependence and radiative 
corrections. It turns out that / s (lGeV) ~ 0.33 GeV for example. In the two-quark scenario 
for /q(980), the decay constants f n , s are related to via 
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f s = f s cos 9, f n = f n sin 6. 
The hard spectator function relevant for B — > f$K decay has the form 
fufs [ l d P ^ . . f 1 d£ f 1 dn 



Hi(f K] 



Fo f °(0)m 2 



B 

for % 



Jo p Jo 4 Jo V 

6, 8. where £ = 1 — ^ and fj 



(2.10) 



,(2.11) 



for i = 1,4, 10 and Hi = for i = 6,8. where £ = 1 — £ and fj = 1 — n. As for the parameters 
agg(K°/o) appearing in Eq. (|2.1|) . they have the same expressions as agg(/oK°) except that the 
penguin function Gk (see Eq. (55) of the second reference in |l6(]) is replaced by Gf n and <& P K by 
For the distribution amplitudes $j , and the annihilation amplitudes, see [22j for details. 
Although the parameters a,(i 7^ 6, 8) and a6,8 r x are formally renormalization scale and 75 scheme 
independent, in practice there exists some residual scale dependence in aj(yu) to finite order. To be 
specific, we shall evaluate the vertex corrections to the decay amplitude at the scale u = nib- In 
contrast, as stressed in 0], the hard spectator and annihilation contributions should be evaluated 
at the hard-collinear scale = \/fIKh with A^ ~ 500 MeV. There is one more serious complication 
about these contributions; that is, while QCD factorization predictions are model independent 
in the m& — > 00 limit, power corrections always involve troublesome endpoint divergences. For 



7 



example, the annihilation amplitude has endpoint divergences even at twist-2 level and the hard 
spectator scattering diagram at twist-3 order is power suppressed and possesses soft and collinear 
divergences arising from the soft spectator quark. Since the treatment of endpoint divergences is 
model dependent, subleading power corrections generally can be studied only in a phenomenological 
way. We shall follow to parametrize the endpoint divergence Xa = Jq(Ix/(1 — x) in the 
annihilation diagram as 

* A = lnte)(l + ^e^) (2.12) 

with pA < 1. Likewise, the endpoint divergence Xjj in the hard spectator contributions can be 
parametrized in a similar way. 



B. Consideration of mixing-induced CP asymmetry 

Consider the mixing-induced CP violation in the decay modes (<p, u, p°, ir°, rf, rj, fo)Ks mediated 
by b —> sqq transitions. Since a common final state is reached only via mixing, hence 

^ _ Qb_ Qk_ A(B° -► MK°) 



p B Pk A(B° -» MK°) 



=r0 



v tt v td V: d V cs A(B - MK 



~ V th V* d VcdV* A(B° - MK°) ■ (2 - 13) 
We shall use this expression for A/ to compute CP asymmetries Sf and Af. For M = V, A(B° — ► 
has the same expression as — A(B° — > t^if ) with the CKM mixing angles X p — ► A* owing 
to CP|F£r°) = -|F#°), while for M = P, A(B° -> PA" ) is obtained from ,4(P° -> PiT°) with 
A p — > A*. If the contributions from V^&V^, terms are neglected, then we will have 

l£ZS ""W ^>""»%=W-» ^S,~W (2,4) 
Note that rjf = 1 for /o-K"s an d 7// = —1 for (0, a;, p , rj', 77, ir°)Ks. 

From Eq. (|2.1j) we see that among the seven modes under consideration, only uKg, p°Ks, tt°Ks 
and rjC^Ks receive tree contributions from the tree diagram b — > sqq (q = u,d). However, since 
the tree contribution is color suppressed, the deviation of —rjfSf from sin 2(5 is expected to be 
small. Nevertheless, the large cancellation between 04 and clq penguin terms in the amplitudes 
of B° — > loK° and B° — ► p°K° render the tree contribution relatively significant. Hence, ASf is 
expected to be largest in the uiKs and p K$ modes. Since the typical values of the effective Wilson 
parameters obtained from Eq. (|2.5|) are 

a 2 « 0.18 - O.lli, 03 ps -0.003 + 0.005i, a 5 ps 0.008 - 0.006i, 

a\ ps -0.03 - 0.02i, a\ ps -0.03 - 0.006i, 

ps — 0.06 — 0.02z, ps -0.06 - 0.004i, (2.15) 

and ps 0.57, it is not difficult to see from Eq. (|2.1[) that 5f lies in the region > 5f > —ir/2 for 
the wi^s and foKg modes, ir > 5f > tt/2 fox and 7r/2 > #f > for the remaining three ones. 

Therefore, based purely on SD contributions, it is expected that ASf > for all the modes except 
for p°K$ and that Af is negative for loKs, foKs and positive for (f)Ks, p°Ks,TT Ks,rj^Ks- 
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C. Numerical results 



To proceed with the numerical calculations, we shall follow 16, 19] for the choices of the relevant 
parameters except for the form factors and CKM matrix elements. For form factors we shall use 
those derived in the covariant light-front quark model and assign a common value of 0.03 for the 
form factor errors, e.g. F Bn (0) = 0.25 ±0.03. For CKM matrix elements, see the unitarity triangle 
analysis in [2^. For definiteness, we shall follow the first reference in [2^ to use the Wolfenstein 
parameters A = 0.801, A = 0.2265, p = 0.189 and fj = 0.358 which correspond to sin 2/3 = 0.723 
and 7 = 63°. We assign 15° error for the unitarity angle 7, recalling that two values 7 = (62l{2)° 



and 7 = (64 ± 18)° are obtained in |25j. For endpoint divergences encountered in hard spectator 
and annihilation contributions we take the default values pa = Ph = 0. We will return back to this 
point below when discussing long-distance rescattering effects. 

The obtained branching ratios for the decays B — > (0, u, p°, ir°, 77', 77, fo)K° are shown in the 
second column of Table [0J while the corresponding CP violation asymmetries Sf and Aj are 
depicted in Table [TTTl In general, our results for branching ratios and direct CP asymmetries are 
in agreement with [l9(. Some differences result from different inputs of the form factors and CKM 
parameters. It is evident that, as far as the central values are concerned, the predicted branching 
ratios by the short-distance (SD) QCD factorization approach are generally too low compared to 
experiment especially for ujK° , p Q K$ and 7r°K° . Note that B(B° — > u>K°) < 10~ 6 is predicted in 



the early QCD factorization calculation 26j. The very large (small) branching ratio for r]'K° (r]K°) 
is understandable as follows. There are two distinct penguin contributions to r]^K°: one couples 
to the d quark content of the rf \ while the other is related to the s quark component of the 7/ ) 
[see also Eq. (|2.1jl ]. If the rj — rj' mixing angle is given by —19.5°, the expressions of the 7/ and 77 
wave functions will become very simple: 

|t/) = —=\uu + dd + 2ss), 1 77) = —=\uu + dd — ss) . (2-16) 
v6 y3 

It is evident that the SD r]K° amplitude vanishes in SU(3) limit, whereas the constructive inter- 
ference between the penguin amplitudes accounts for the large rate of rj'TC . In reality the 77 — 77' 



mixing angle is —(15.4 ± 1.0)° |2]J], but this does not affect the above physical picture. 



Owing to the large cancellation between the 04 and a§ penguin terms, the main contribution 
to the decay S° — > fo~K° arises from the penguin diagram involving the strange quark content of 
/o(980), namely, the term with the scalar decay constant f s . Consequently, the maximal branching 
< 10" 
Tab! 

by BaBar with the result 



ratio 9.9 x 10 e occurs near the zero mixing angle. The result of B(B° — > foK°) = 8.1 x 10 e 
shown in Table ITT1 corresponds to 9 = 150°. Note that the decay i?° — ► fo(980)K° was measured 



B(B° -> / (980)K° -» tt + tt-K°) = (6.0 ± 0.9 ± 1.3) x 10" 6 . (2.17) 
The absolute branching ratios for B — > foK depends critically on the branching fraction of 
/o(980) — > 7T7T. We use the results from the most recent analysis of [3| to obtain £>(/o(980) — * 
vrvr) = 0.80±0.14 and B(B° -» / (980)i^°) = (11.3±3.6) x 10~ 6 as shown in TableHH 1 In short, al- 



1 For comparison, the world average of the branching ratio for B — > foK — > tt + tt K is (8.49^11) x 10 6 
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TABLE II: Short-distance (SD) and long-distance (LD) contributions to the branching ratios (in 
units of 10 -6 ) for various penguin-dominated modes. The first and second theoretical errors cor- 
respond to the SD and LD ones, respectively (see the text for details). The world averages of 
experimental measurements are taken from |l7j |. 

SD SD+LD Expt 



B — > loK 2.0^^'3 5.6^^*2^2*1 - 1 - u : ' 

B°^p K° 2.8±H 5.2±? : i±f : | 5.1 ±1.6 

B^^iTf 42.1^ 69.4±ii;|±f§;| 68.6 ±4.2 



jf^ V K° i-s±S:i i-s±Si±8:S <2.o 

B°^7T K° 5.8±l;f 9-6±i:g±l:g 11.5 ±1.0 

B° ^ f K° 8.1±|;^ 8.1±l;^±g:g ** 11.3 ±3.6 



a Only the intermediate states if p + and 7r + if * are taken into account; see Sec. Ill for details. 

though the predicted branching ratios of (<j), u, p°,rf, tt )K° are consistent with the data within the 
theoretical and experimental uncertainties, there are sizable discrepancies between the SD theory 
and experiment for the central values of their branching ratios. This may call for the consideration 
of subleading power corrections such as the annihilation topology and/or FSI effects. 

In Tables ITT1 and ITTT1 we have included the SD theoretical uncertainties arising from the variation 
of the unitarity angle 7 = (63 ± 15)°, the renormalization scale p from 2rrib to m&/2, quark masses 
(especially the strange quark mass which is taken to be m s {2 GeV) = 90 ±20 MeV) and form factors 
as mentioned before. To obtain the SD errors shown in Tables II and III, we first scan randomly 
the points in the allowed ranges of the above four parameters in two separated groups: the first 
one and the last three ones, and then add each error in quadrature. For example, for the decay 
B° - n'K s we obtain 2B = (42.1^^) x lO" 6 , A = 1.77l°;?°l°Ju% and S = 0.737t8»jj& 
where the first error is due to the variation of 7 and the second error comes from the uncertainties 
in the renormalization scale, the strange quark mass and the form factors. 

From Table |TTT| we obtain the differences between the CP asymmetry Sj D induced at short 
distances and the measured Sjj^ Ks to be 

A cSD _ n n9 +0.00 a q SD _ n 1 9 +0.05 a qSD _ n nn+0.03 

A qSD _ n nR+ 02 a ? SD — n m + - 00 a <? 5D - n n7+ - 02 

^^it°K s — U - UD -0.04> ^^v'Ks — U ' Ui -0.04; ^^-qKs ~ U - U '-0.04> 

AS^ Ks = 0.021°;°°, (2-18) 

where the experimental error of Stmk s 1S n °t included. Our results for AS^ S and AS r 5^- s are 
consistent with that obtained in [19j. 2 As expected before, the uKg and p°K$ modes have the 



Q and hence B(B~ -» f K~) « (15.911;?) x 1(T 6 . 
2 Note that unlike [19| we did not include the theoretical uncertainties arising from power corrections. 
Otherwise, there will be a double counting problem when considering LD rescattering effects. 
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TABLE III: Short-distance (SD) and long-distance (LD) contributions to the time-dependent 
CP asymmetry. The first and second theoretical errors correspond to the SD and LD ones, respec- 
tively (see the text for details). 



Final State 




~ n f S f 






A f (%) 






SD 


SD+LD 


Expt 


SD 


SD+LD 


Expt 



<f>K s 

loK s 

P°K S 

V'K S 

j]K s 

TT°K S 

f K s 



0.747; 



-0.002 
-0.039 

o.85oi8:°g 

n fi orr+0.028 

n 7 o 7 +o.oo2 

u - '°' -0.038 
-0.017 
-0.044 
n 7CV+0.018 
U -' 8 ' -0.044 

u - '* y -0.039 



0.759; 



0.793 



-0. 

-0. 

0.73618; 

0.76118; 
0. 

-0. 
0. 
0. 
+0. 



0.725 
0.802 
0.770 



0.74918 



007+0. 

041- 0. 
022+0. 

035- 0. 
071+0. 
079-0. 
004+0. 

036- 0. 
025+0. 
046-0. 
006+0. 

042- 0. 
002+0. 
039-0. 



005 
006 
025 
014 
073 
100 
005 
003 
002 
004 
015 
019 
a 




0.35 ±0.20 

^+°- 30 
u ' oo -0.32 



0.43 ±0.11 



0.34 



+0.27 
-0.29 



0.39 ±0.26 0.77 



1 4+0-3 
L - q: -0.5 
7 q+3.5 
" ' - -2.6 

9.0111 

j+0.4 
'-0.4 
+5.1 
-2.0 
+2.1 
-1.1 
+0.13 
-0.10 



9 fi +0.8+l.l 
" z -°-1.0-0.9 



-13-2111 
46.6±H;fl 



1.1 
-6.1 
-3.4 



2.1 



+0.6- 
-0.3- 



2.1 
;-2.6 
-10.8 
-5.9 
0.5 
0.2 



o 7+4.4+1.4 
°- '-1.8-2.4 
o 7+I.8+2.O 
°- '-2.0-0.4 
+0.14+0.01 
-0.09-0.01 



0.75 



4± 17 
48 ±25 

4±8 

-8± 14 
-14 ±22 



"Only the intermediate states K p + and ir + K* are taken into account (see Sec. Ill for details). This 
means that the prediction of the LD effects on the JqKs mode is less certain. 



largest deviation of Sf from the naive expectation owing to the large tree pollution. In contrast, 
tree pollution in rj'Ks is diluted by the prominent ss content of the if. As for direct CP violation, 
sizable direct CP asymmetries are predicted for uK$ and p°K$ based on SD contributions. 



III. LONG-DISTANCE CONTRIBUTIONS 



As noticed in passing, the predicted branching ratios for the decays ~B° — > (4>,u>, p° ,ir°)Ks by 
the short-distance QCD factorization approach are generally too low by a factor of 2 compared to 



experiment. Just like the pQCD approach |29J where the annihilation topology plays an essential 
role for producing sizable strong phases and for explaining the penguin-dominated VP modes, 
it has been suggested in that a favorable scenario (denoted as S4) for accommodating the 
observed penguin-dominated B — > PV decays and the measured sign of direct CP asymmetry in 
B° — » K~ir + is to have a large annihilation contribution by choosing pa = 1, 4>a = —55° for PP, 
(j)A = —20° for PV and cpA = —70° for VP modes. The sign of 4>a is chosen so that the direct 
CP violation Aj^- n + agrees with the data. However, there are at least three difficulties with this 
scenario. First, the origin of these phases is unknown and their signs are not predicted. Second, 
since both annihilation and hard spectator scattering encounter endpoint divergences, there is no 
reason that soft gluon effects will only modify pa but not pn- Third, the annihilation topologies 
do not help enhance the 7r vr and p°7r° modes; both pQCD and QCDF approaches fail to describe 
these two color-suppressed tree-dominated modes. As stressed in [l9(, one would wish to have an 
explanation of the data without invoking weak annihilation. 

As shown in ^jj, final-state rescattering can have significant effects on decay rates and CP vi- 
olation. For example, the branching ratios of the penguin-dominated decay 4>K* can be enhanced 
from ~ 5 x 10 -6 predicted by QCDF to the level of 1 x 10~ 5 by FSIs via rescattering of charm 
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intermediate states jlil ]. Indeed, it has been long advocated that charming-penguin long-distance 
contributions increase significantly the B — > Ktt rates and yield better agreement with experiment 



30 



31]. The color-suppressed modes D°tt°, 7r°7r° and p ir° in B decays can also be easily enhanced 



by rescattering effects. Moreover, large nonperturbative strong phases can be generated from the 
final-state interactions through the absorptive part of rescattering amplitudes. We have shown ex- 
plicitly in that direct CP-violating partial rate asymmetries in K~ir + , p + -K~ and tt + tt~ modes 
are significantly affected by final-state rescattering and their signs, which are different from what 
is expected from the short-distance QCDF approach, are correctly predicted. In order to avoid 
the double-counting problem, we will turn off the LD effects induced from the power corrections 
due to non- vanishing pa and pn; that is, we set pa = pu = and 4>a = 4>h = 0, therefore it is 
important to note that we are not adding FSI on top of QCDF. We wish to stress that in principle 
LD rescattering effects can be included in the framework of QCDF, but that requires modelling of 
power corrections and, in particular, one may then need to adopt non- vanishing values 
of pa, Ph, 4>a and 4>h 0|, as mentioned above. In this work, we are providing a specific model for 
final-state rescattering to complement QCDF. 

Besides direct CP violation, the mixing-induced CP asymmetry Sf also could be affected by 
final-state rescattering from some intermediate states. When the intermediate states are charmless, 
the relevant CKM matrix element is V u bV* s ~ ^4A 4 i?& e ~* 7 which carries the weak phase 7. In 
general, the charmless intermediate states will essentially not affect the decay rates but may have 
potentially sizable effect on Sf, whereas the charm intermediate states will affect both the branching 
ratios and Sf. 



A. Final-state rescattering 

At the quark level, final-state rescattering can occur through quark exchange and quark annihi- 
lation. In practice, it is extremely difficult to reliably calculate the FSI effects, but it may become 
amenable to estimate these effects at the hadron level where FSIs manifest as the rescattering 
processes with s-channel resonances and one particle exchange in the t-channel. In contrast to 
D decays, the s-channel resonant FSIs in B decays is expected to be suppressed relative to the 
rescattering effect arising from quark exchange owing to the lack of the existence of resonances at 
energies close to the B meson mass. Therefore, we will model FSIs as rescattering processes of 
some intermediate two-body states with one particle exchange in the i-channel and compute the 
absorptive part of the rescattering amplitude via the optical theorem 0] . 

Given the weak Hamiltonian in the form Hyj = J2i^iQi> where Aj is the combination of the 
quark mixing matrix elements and Qi is a T-even local operator (T: time reversal), the absorptive 
part of final-state rescattering can be obtained by using the optical theorem and time-reversal 
invariant weak decay operator Qi. From the time reversal invariance of Q (= U T Q*U ] T ), it follows 
that 

(i; out\Q\B; in)* = £ S^j; out|Q|B; in), (3.1) 
j 

where Sij = (i;out|j;in) is the strong interaction S-matrix element, and we have used 
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f/rlout (in))* = |in(out)) to fix the phase convention. Eq. (|3.1|) implies an identity related to 
the optical theorem. Noting that S = 1 + iT, we find 

2Abs{i;out\Q\B;m) = ^ T^{j; out\Q\B; in), (3.2) 



where use of the unitarity of the 5-matrix has been made. Specifically, for two-body B decays, we 
have 



1 



Abs M(p B - p a pb) = - J2[ n l=i 



(2^2E k 



x S\ Pa +p b - Qk)M( PB - fe})T*( M -> {%}). (3.3) 
fe=l 

Thus the optical theorem relates the absorptive part of the two-body decay amplitude to the sum 
over all possible B decay final states {qk}, followed by the strong rescattering {q/.} — ► p a Pb- In 
principle, the dispersive part of the rescattering amplitude can be obtained from the absorptive 
part via the dispersion relation 

V f°°AbsA(s') 



o. r f°° AbsA(s') , , 
VisA{m B ) = — / — VW. 

7T Js s — 



(3.4) 



Unlike the absorptive part, it is known that the dispersive contribution suffers from the large 
uncertainties due to some possible subtractions and the complication from integrations. For this 
reason, we will assume the dominance of the rescattering amplitude by the absorptive part and 
ignore the dispersive part in the present work. 

The relevant Lagrangian for final state strong interactions is given by 



where 



--Tr[F^(V)F^(V)]+ig V ppTr(V^P 8 „P) + gvvpe^Ti{d^V u d a VpP) 



(3.5) 



1 



4 

+igvpp 



l K- d + p + »(K° d „K- - V2tt° 8^-) + K*~ 



+K*-(^v' q d»K+ + K+ d X) + K *°(^' Q d„K° + K° 5^) + . . . 



+9WP e 



with P and V b eirrg the usiicil pseudosc&lcir and vector niultiplets, respectively, 
igvWm V u }/2, ffj, u = dfj,V u - duV/j, and 



(3.6) 



1 



C h = - igD. D p(p i d>>PijDf -D^d^PijD^-^gD^pe^pD^P^ d a D*f 
~ ig DDV D\ \ Di(V»ti - 2f D , DV e^{d^)){D\ d a D*^ - Df d a D>) 



+ ig D *D*vB^ <9 M D*J>(V»)) + Aif D , D , v D*^V u - d v V%D*J. 



(3.7) 
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K°( Pi ) 



D^- (pi 




4>ipz) 



D { * ]± {k) 



K°( P4 



(b) 



FIG. 1: Final-state rescattering contributions to the B° — > (pK° decay. 



Only those terms relevant for later purposes are shown in L\ and the convention g 0123 = 1 has been 
adopted. For the coupling constants, we take q p kk — Qp-mx — 4. 28, 3 g^KK = 4.48, y/2 gyyp = 
16 GeV -1 In the chiral and heavy quark limits, we have [3^ ] 



9d*d* 



9D*Dtt 



9ddv = 9d*d*v 



P9V 



Id*dv 



Id*d*v ^9v 



(3.8) 



■rr ) JU-UV — it: 1 

m D * V2 rn D * a/2 

with f n = 132 MeV. The parameters gy, f3 and A (not to be confused with the Wolfenstein 
parameter A) thus enter into the effective chiral Lagrangian describing the interactions of heavy 
mesons with low momentu m lig ht vector mesons (see e.g. [3^] ). The parameter gy respects the 
relation gy = m p / = 5.8 |23j. We shall follow |3j| to use (3 = 0.9 and A = 0.56 GeV -1 . The 
coupling gr>*D-K has been extracted by CLEO to be 17.9 ± 0.3 ± 1.9 from the measured D* + width 



B. B° — > 4>Ks as an example 

We next proceed to study long-distance rescattering contributions to the b — > sqq transition- 
induced decays B° — > (c/>, u>, p°, rf\ 7r°, Jq)Ks- To illustrate the calculations of rescattering ampli- 
tudes , we shall take the 4>K$ mode as an example. Its major final-state rescattering diagrams are 
depicted in Fig. 1. 

The absorptive parts of the B° — > K~ p + — > (j)K° amplitude via the K^, K*^ exchanges are 
given by 

,,,„-+ r,±, 1 f d 3 pi d 3 p 2 . A(B°^K-p+) 
Abs{K p^;K^) = - / - {^) o^[p B - Pi- P2) 



2 J (2tt) 3 2E 1 (2ir) 3 2E 2 v ' ^ 1 * ' 2e\ ■ p B 
^2 2e 2 • PB (-1i)94>KK{£t ■ Pi) —. ^ t,m 2 K ^ 2i 9pKK {£2 ■ Pi) 



x , _ . . • ■ ,., .> . ... - , 



9* v 1 {^lE^^lI^io \W( A(B°^K-p^ 

2 ^3 • PB X - / jr-r^— ,„ (2tt) S ( PB ~ Pi ~ p 2 ) 



2 J (2tt) 3 2Ei [2tt) 3 2E 2 k ' ^ 2e* • p B 

F 2 (t,m K ) m (i) / (jpx -p 2 )(P2 -PA) 



xig+KK /' 2 g P KK{A\ LI - A£>) -pi ■ pa + ^' y 2 . (3-9) 



t — m 2 K 



The pirir coupling defined here differs from that in 0] by a factor of l/\/2. 
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(2tt) 3 2E 1 (2tt) 3 2E 2 v ' ^° ^ 2e\ ■ p B 

(k^k v ^ 
-9^ H — 

F 2 (t,m K *) 



A' * 



2e 3-PBX- /n . Uop /0 .,,„ P (2vr) 5 (p B - Pl - pa) 



2 7 (2^)32^1 (2vr) 3 2^ 2 2e^-p B 
xg<f,K*K -, (t,m f^ g pA ^^(-24 2) m|), (3.10) 

where the dependence of the polarization vector in the amplitude A(B° — » K~ p + ) has been ex- 
tracted and k = pi — p$. In order to avoid using too many dummy indices, we have defined 
[A,£,C,Z>] = t a ^A a B^CPD\ [A,B,aj.] = e aP ^A a B^C< and so on for later convenience. 
Moreover, we have applied the identities |l5( 

p lfl = P^ ] +q^\ 

PlpPlu = g^A® + P^Af + (P„q u + q^Af + ^^4 2) > (3-11) 

under the the integration 

d 3 Pi d*p 2 ftt _^. 



(2^)32^ (2^)32^2 ^ ^ (PS ~ Pl ~ P2)f{t) X (3 ' 12) 
with P = pz +pa, q = P3 — pa- These identities follow from the fact that the above integration can 
be expressed only in terms of the external momenta P3,P4 with suitable Lorentz and permutation 
structures. The explicit expressions of Af = Af(t,m B ,mlmimimj) can be found in Q. 

Before proceeding it should be stressed that we have applied the hidden gauge symmetry La- 
grangian Eq. ()3.6|) for light vector mesons and the chiral Lagrangian Eq. Q3.7J1 . based on heavy 
quark effective theory (HQET) and chiral symmetry, for heavy mesons to determine the strong 
vertices in Fig. 1. This requires that the involved light pseudoscalar or vector mesons be soft. 
However, the final state particles are necessarily hard and the particle exchanged in the t channel 
can be far off shell, especially for the t-exchanged D meson. This is beyond the applicability of the 
aforementioned chiral perturbation theory and HQET. Therefore, as stressed in 15|, it is necessary 
to introduce the form factor F(t, m) appearing in Eqs. ()3.9|) and (|3.1Uj) to take care of the off-shell 
effect of the t-channel exchanged particle and the hardness of the final particles. Indeed, if the 
off-shell effect is not considered, the long-distance rescattering contributions will become so large 
that perturbation theory is no longer trustworthy. For example, since B — > D S D is CKM doubly 
enhanced relative to B — * Kit, the rescattering process B — » D S D —* Kir will overwhelm the initial 
B — > Kir amplitude. Hence, form factors or cutoffs must be introduced to the strong vertices to 
render the calculation meaningful in perturbation theory. 
The form factor F(t, m) is usually parametrized as 

( A? - m 2 \ n 

F(t,m)= A2 _ t , (3.13) 
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normalized to unity at t = m 2 with A being a cutoff parameter which should be not far from the 
physical mass of the exchanged particle. To be specific, we write 



A 



+ V^QCD, 



(3.14) 



where the parameter r\ is expected to be of order unity and it depends not only on the exchanged 
particle but also on the external particles involved in the strong-interaction vertex. The monopole 
behavior of the form factor (i.e. n = 1) is preferred as it is consistent with the QCD sum rule 
expectation [33] . Although the strong couplings are large in the magnitude, the rescattering ampli- 
tude is suppressed by a factor of F 2 (t) ~ m 2 AQ CD /i 2 . Consequently, the off-shell effect will render 
the perturbative calculation meaningful. Moreover, since in the heavy quark limit t ~ m B , the 
final-state rescattering amplitude does vanish in the ms 



Likewise, the absorptive part of the B 
given by 



K* 



7T 



— > 00 limit, as it should. 
(j)K° amplitude via the K*^ exchange is 



Abs (K'-jr+iK**) 



(27r) 3 2£a (27r) 3 2£ 2 
9V 



(2tt) 5 (p B -P1-P2 



A(B° 



p B i 



4if 



2 



in. 



K* 



l^ei-ps) 

F 2 (t,m K *] 
t 



2eJ • PB 

£l • 4(Pl + P3)ti + £l M ( 2 Pl • £ 3)] 



m 



K* 



2e* 3 ■ p B x 



1 f d 3 pi d 3 p 2 



2 J (2vr) 3 2^i (2vr) 3 2£ 2 



-i)gK*Kir(P2 +PA) 



(27r) 5 (p B ~P\-P2 



A{B' 



K 



(3.15) 

-7T + ) 



2e\ ■ p B 



x 9v 



F 2 (t,m K * 
2(t 



m 



■ 9K*Ktt < 



K* 



A^ + A 



(i) 



(A 



(i) 



A 



2 



in 



X [P2-P3 



(pi -P2){P2 -P3] 



+ 



(i) 



A 



(i) 



1 + (A 



(i) 



K* 



A^ ] ) Pl T 2 



(Pi +P3) ■ (P2 +Pi) + 



mi 



in 



m 2 )(m 



m A 



K" 



2(4 1) -^) 



0-h 



m 



P2- 



K* 



niA 



in 



+ Pi- 



nt 



K* 



m\ 



m% 



K* 



in 



K* 



P2 -Pl- 



Pl -P2 



Note that since the ttKK vertex is absent, there is no contribution from the K + exchanged particles. 

The B° — > (pK° decay also receives contributions from charmless VV modes. The leading candi- 
date is the K*~ p + mode via the p-wave configuration. However, as we have checked numerically, its 
amplitude is one or two orders of magnitude smaller than those from the previous two charmless VP 
modes and its effect on Sf is quite small. Given this, we will not go into any further detail on the 
rescattering from charmless VV modes. Thus far we have only considered the contributions where 
the two intermediate states originating from the weak vertex are on shell. There are additional 
contributions where one of the mesons coming from the weak vertex and the exchanged meson in 
the t-channel are on shell. For example, in the diagram of Fig. 1(a), we can set K~ and K + on shell 
while keeping p + off shell. This corresponds to the 3-body weak decay B° — » K + K~K° followed 

' " (j)K° where K° behaves as a spectator. However, there are 



by the strong rescattering K + K K 
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-ftO 



many possible pole contributions to B 
inferred from Fig. 1(a), one can also have B — > B S K 



K + K~K°. In addition to B° 



K-p 



K + K~K° 



as 



K + K~ K° , for example. In the present 
work, we will only focus on two-body intermediate state contributions to the absorptive part. Since 
the analogous 3-body contributions do not occur in Fig. 1(b) and since Fig. 1(a) is CKM doubly 
suppressed relative to Fig. 1(b), it is safe to neglect the additional three-body contributions for 
our purposes. 

We next turn to the FSI contribution arising from the intermediate states D^*\ Note that 
only the p-wave configurations of the systems can rescatter into the (f>K° final state. The 

absorptive parts of B° — ► D*~D + — > 4>K° amplitudes via D* exchanges, B — > D*~D + , D*~D* + — > 
(f)K° amplitudes via D s and D* exchanges, are given by 

d 3 Pi d 3 p 2 /0 _, 4r4 ,„_ ,A(B°^D*-D+) 



•Abs (D*~D + ;D* ± ) = 1 / ^ ^ (2,)^ - Pl - p 2 ): 



Ate (D7I>* + ; Of) 



Abs (D~D* + ;D* ± ) 



X^dpuPl ■ 4 + P3p£*3n ~ 4 P P3v) 

d 3 pi d 3 p 2 



~i 9d*dk 



k^k u 



m 



Di 



2e\-P 



Plv 



= H-PB x\j 



(27r) 4 (5 4 (p jB - Pi -P2) 



(2tt) 3 2£i (27r) 3 2£ 2 



A(B° 



2e\-P 



D*-D+) F*{t,m D .) 

— ^iDtD*^— 9 9D*DK 



m 



D: 



Pl -P2 



-Pi 



,k-p 4 

VA - 



m 



+ 



i - ( 4 1} 

1 - (A? 



A 



(l)xfc Vg4 



771 



P2 "P3 



Pl • P2 Pl • P3 



mt 



<i 3 pi d 3 p*2 



P3 • P4 ~ 



P3 • k k -p4 



rn 



D< 



A 2 



2£3 ' PS x 2 



(2vr) 5 (p B - pi -p 2 ) 
i — m 



2e* • P 



(2vr) 3 2£i (2tt) 3 2E 2 
}_^2e 2 -P 2ig DsDs4> e 3 ■ pi — — — 3 — {-i)9d*d s k Pa • £2 



(27r) 4 (5 4 (p jB - Pi -p 2 ) 



A(S° 



D~D* 



(2vr) 3 2£i (2ir) 3 2E 2 

L 1 — m Dg \ 

1 [ d?Pi d 3 p 2 44 

2 7 (2^2^ (^2^ (2 " } 5 " Pl " ^ 



+ 



2e\-P 

Pl • P2 P2 • P4 



A(B° £>7l>* + ) 



x^2 £ *-P4z/ D * Ds 

A 2 



F 2 (t,m D * 



m 



-WD*D*K 



Di 
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x[P3,£3,Pi, M ]( -g^ + ^- 



m 



Dt 



2e 3 ' Pb x 2 

x (-8) fD*D s <f> 



1 /" d 3 pi d 3 p 2 



2 7 (2vr) 3 2 J Ei (2tt) 3 2 j E 2 



(27T) (5 (p B -P\-P2) 



A(B° -► DjD*+) 
2eS • P 



t — m 



, 9d*d%k mlAf\ 

D* 



Abs (D*~D* + ;Df) = \ J ( ^ ( ^ (2,)^ - Pl - p 2 )(^, „ P,p 2 ]) 

F 2 (t,m Da 



y: ere? {tWb-d. 



Ai,A2 



{-i)9D*D 3 K [P3,£3'Pl' £ lb4 • £2 



H'Pbx^J 



1 /■ <i 3 pi d 3 p 2 



(27r) 3 2£i (2tt) 3 2£ 2 
x(-4ic)/ g * Ds0 i? / t,m 2 P ^ g D * Ds ^m^{ 2) ) 



(27r) 4 ,5 4 (p B -p 1 -p 2 ) 



t — m 



Abs (D*~D* + ; D*^) = 



1 /" (i 3 px d 3 p 2 



(2tt) (5 (p B -pi -p2)(«c[ 5 , /3 ,P,p 2 ]) 



2 7 (27r) 3 2 J Bi (2tt) 3 2.E 2 

*<5 */3 / , -x r F 2 (t,m D *) . 

X £ 1 £ 2 (-tyfD'DSt-j— — 2^WD*D* 3 K 

Ai,A 2 1 

xeKcrgp^pi ■ £3+ P3p£3 M - £ 3pP3 M ) I -ff^H — ) [£2,P4,P2, v ] 

s 

If d 3 pi d 3 p 2 
2^3 ■ Pb 2 J 



x(-4ic)/n» D » 



(2tt) 3 2£i (2tt) 3 2£ 2 
F 2 (t,mn*) 



(27t) 4 ,5 4 (pb-Pi-P2) 



i-m 2 



-9d*d*k* 



x - 4 1} )(Pi • - pi • p 2 p 2 • p 4 ) + m^ ; ], 



2 /I (2)l 



(3.16) 



where the dependence of the polarization vectors in A(B° — »■ D* D + , D s D* + ) has been extracted 
out explicitly and a = qd*d*v /(2/d*d*v) and the i? — ► D*~D* + decay amplitude has been denoted 



as 



A(5 D* s ( Pl , X 1 )D*( P2 , A 2 )) = efe^ (a^ + + *c [ M , „, P,p 2 ])- 



(3.17) 



In order to perform a numerical study of the above analytic results, we need to specify the short- 
distance A(B — > £>s amplitudes. In the factorization approach, we have 



V2 



V A VZa 1 f D .m D .F* D (m 2 D .)(2e* D . ■ PB ), 
V cb V c * sai f Ds m D ^ D \m 2 Ds )(2s*^ -p B ), 



.Gf , 



A(B ^ D*D*) S D = -i^V a V£aif D *m D z(mB + m D *)e*g,£D, 



Af D * (m 2 D * )g lw 



2A^{m Dt ) 2V^'(m 2 m ) 



(m B + m D *) 



^PB^PBu ~ I 



{tub + mo*) 



2 e [tvapPBPD* 



(3.18) 
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Since the phase of the parameter a\ originating from vertex corrections [see Eq. (|2.5|) ] is very small, 
one can neglect the strong phase of the short-distance amplitudes. 

The long-distance contribution to the B° — > ujK decay can be preformed similarly. Due to 
the absence of the PPP vertex and the G-parity argument, the number of FSI diagrams from 
charmless intermediate states is greatly reduced compared to the previous case. For example, the 
K~ p + intermediate state does not contribute to the K°u amplitude (through t-channel it and 
p exchanges) as both KKtt and ppuj vertices are forbidden. In fact, there is only one relevant 
rescattering diagram, namely, FT K*~tt+ K uj via p exchange, arising from charmless 
intermediate states and the corresponding absorptive part is given by 

Abs(K ir^;p ) = - / f .^.„ frt „ n „ (2tt) 5 (p B ~ Pi ~ p 2 ) 



2 J (2tt) 3 2E 1 (2tt)32£ 2 K ' K1 D 1L 2e\-p B 
x ^2e\ -p B (ig P K*K)(iy/29u;(m) bl> £ l> fc >Mp>f>.P4)4] 

\ m K * J t — m A p 

= 2£4 • PB x 2 J WWm WWe 2 (27r) 6 iPB - pi ~ P2) — WV B — 

x (-2^2) g pK * K F (t ' mp) g^mlAf . (3.19) 
(t — nip) 

Furthermore, we have checked numerically that the B° — ► K*~ p + — > K°lo contribution is small. 
Hence we will skip further detail on the rescattering from charmless VV modes. 

The rescattering amplitudes of B° -> D*- D + , D~ D* + , D*-D* + -> K°u; amplitudes via D, D* 
exchanges can be evaluated in a similar way. The analytic expressions of their absorptive parts are 
similar to those for cpK°. For example, 

Abs(DrD+-D±) = 2SI-PBX 1 - I (2nm PB - Pl - P2 ) A(]f ~* D7 °* +) 

v s ; 4 F 2 J (2vr) 3 2.Ei (2n) 3 2E 2 2e* 2 ■ P 

x |V2g D|DX f _ m ^ 9dDuj {l-A\>-A v 2 ')^-p 2 -p 3 + -5 J j, 

(3.20) 

is the same as Abs{D~ D* + ] Df) in 1)3.16)1 after the replacements qd*d s k — > 9d*dk, fD s D s <t> — > 
Iddw/ V^2, Pi P2, P3 ^ P4 5 ^1 — ^ 1 — ^1 — ^2 and a suitable replacement of the source 
amplitude (note that the replacement of momentum should not be performed in Aj). Likewise, we 
can obtain Abs(D*~D+;D* ± ), Abs(D~ D*+; D*±), Abs(D*~ D*+ , dW±) from Abs(DjD*+; D*±), 
Abs{D*-D + ;D* ± ), Abs(D*-D* + , Di* ]± ) in (l3~TB|) . respectively, with ^j 2) being unchanged and an 
additional overall minus sign for D*~D* + contributions. This similarity is by no means accidental; 
it follows from the so-called CPS symmetry, i.e. CP plus s «-> d switch symmetry. A similar but 
more detailed discussion is given in 0] for the case of B — > <pK * decay. 

The final-state rescattering contributions to other penguin-dominated decays B° — ► 
(p° ,rj('\ fo)Ks can be worked out in a similar manner. The dominant intermediates states for 
each decay modes are summarized in Table ITVl 
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TABLE IV: Dominant intermediate states contributing to various final states. For the rjKg final 
state, the intermediate states are the same as that for rj'Ks- 



Final state 


Intermediate state (exchanged particle) 


4>K S 


K*-ir + (K*), K~p + (K,K*), 




D* S -D+(D* S ), D-D*+(D S ,D* S ), D*~ D* + (D S , D*) 


loK s 


K*-ir+(p), 




D*-D + (D,D*), D-D* + (D*), D*~D* + (D,D*) 


p°K s 


K*-7T + (p), 




D*-D+(D,D*), D-D*+(D*), D*-D* + (D,D*) 


rfK s 


r]'K°(K* ), K*-p+(K,K*) 




D-D+(D*,D* S ), D*-D*+(D, D s , D*,D* S ) 


ir°K s 


rj'K°(K* ) 




D~D+(D*), D*-D* + (D,D*) 


f K s 


p+K-(K,p), tt + K*-(tt) 




D*-D+(D,D*), DjD*+(D*) 



C. Results and discussions 

Writing A = A SD +iAbsA LD with AbsA LD obtained above and the form factors given in the 
results of the final-state rescattering effects on decay rates, direct and mixing- induced CP violation 
parameters are shown in Tables [H] and II I II As pointed out in > the long-distance rescattering 
effects are sensitive to the cutoff parameter A appearing in Eq. (|3.14f) or 77 in Eg. (|3.14|) , Since we do 
not have first-principles calculations of r], we will determine it from the measured branching ratios 
and then use it to predict the CP- violating parameters Af and Sf. As shown in 15], r] = 0.69 for 
the exchanged particles D and D* is obtained from fitting to the B — > Kir rates. We take r/ = 0.85 
for the exchanged particles and to fit the data of B° — > rj'lC rates and 77 = 1 for other 
light exchanged particles such as ir,K,K*. As for the VP modes, namely, 4>K$, u)K$ and p°Ks, 
we take i]D (s ) = i]D* = 0.95, 1.4 and 1.5 respectively in order to accommodate their rates 4 . Note 
that the result ilD {s) = Vd* ^ = 1-5 for p K$ is consistent with the value of r\r> = f]D* = 1-6 obtained 
for B — > pn decays within SU(3) symmetry. For the decay if — > foKs, since only the strong 
couplings of /o to KK and ttit are available experimentally, we shall only consider the intermediate 
states K~ p + and -k + K*~ . The LD theoretical uncertainties shown in Tables ITTl and ITTTl originate 



4 Note that a monopole momentum dependence is used for the form factors throughout this paper [see 
Eq. l|3.13fl ]. If a dipole form of the momentum dependence is used for D*-exchange diagrams as in 
Q, we will obtain B SD +™(<f>K°) = (6.1±?;g±g;l) x 1CT 6 , Sfg LD = 0.723^1°^, ^ LD = 
-0.070™ 5 7 ± 0.019, B^°{u>K°) = (2.5+?:™) x 10~ 6 , S*£ LD = 0.847±g:S5±S:gS and A S J+ LD = 

n 01 q+0. 034+0. 027 
U.UlO_ 078-0 .029- 
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from three sources: an assignment of 15% error in Aqcd, the measured error in the coupling 



9d* Dn = 17.9 ± 0.3 ± 1.9 35] and 5% error in the form factors for B to transitions. They 



are obtained by scanning randomly the points in the allowed ranges of the above-mentioned three 

parameters. The calculations of hadronic diagrams for FSIs also involve many other theoretical 

uncertainties some of which are already discussed in [15j . Prom Tables [H] and IIIII it is clear that 

the SD errors are in general not significantly affected by FSI effects and that LD uncertainties are 

in general comparable to the SD ones. 

We see from Table [HI that final-state rescattering will enhance the decay rates of 

4>Ks, ojKsi p°Ks, rj'Ks and ir°Ks but it does not affect the rjKs rate. The seemingly large disparity 

between rj' Kg and rjKs for FSIs can be understood as follows. There are two types of exchanged 

particles in the rescattering processes, namely, and Ds (see Table HVj) . The former (latter) 

couple to the d (s) quark component of the r/') . Since the 7/ and rj wave functions are approximately 

given by Eq. ()2.16j) . it is clear that the rescattering amplitudes due to the exchanged particles D^*> 
(*) 

and Ds interfere constructively for the r] Kg production but compensate largely for r/Kg- 

It should be stressed that although we have used the measured branching ratio of r/'K to fix the 
LD contributions and the unknown cutoff parameter rj, there exist some other possible mechanisms 
that can help explain its large rate. For example, the QCD anomaly effect manifested in the two 
gluon coupling with the rf may provide a dynamical enhancement of the r/'K production j^]. 
And it is likely that both final-state rescattering and the gluon anomaly are needed to account for 
the unexpectedly large branching fraction of rj'K. Note that both contributions carry negligible 
CP-odd phase. Hence, whether the anomalously large branching ratio of rj'K comes from the 
QCD anomaly and/or from final-state rescattering, it will be very effective in diluting the uu tree 
contributions and rendering AS v /k s small. 

We are not able to estimate the long-distance rescattering contributions to the foK$ rate from 
intermediate charm states due to the absence of information on foDD and JqD^D^ couplings. 

Since the modes tvKs, p°Ks,4>Ks and tt°Ks receive significant final-state rescattering contri- 
butions, it is natural to expect that their direct CP asymmetries will be affected accordingly. It 
is clear from Table UTTI that the signs of Af in the last two channels are flipped by final-state in- 
teractions. As for the mixing-induced CP violation Sf, we see from Table HTH that loK s and p°K s 
receive the largest corrections from final-state rescattering, while the long-distance correction to 
<j)Ks is not as large as what was originally expected. The underlying reason is as follows. The 
mixing CP asymmetries S^k s an d & or s deviate from sin 2(3 as they receive contributions from 
the tree amplitude. The contribution from the tree amplitude is relatively enhanced as the QCD 
penguin amplitude is suppressed by a cancellation between penguin terms {\a^ — r x ae\ <C |a,4 1 , \ae\). 
Final-state rescattering form D^D^ states with vanishing weak phases will dilute the SD tree 
contribution and bring the asymmetry closer to sin 2/3. For the (j)Ks mode, the asymmetry 
is slightly greater than sin 2/3. With rescattering from either charmless intermediate states, such 
as K*ir, or from Z)W states, the asymmetry is reduced to S^Ks — 0.73. When both charmless 
and charmful intermediate states are considered, the asymmetry is enhanced to S^k s — 0-76 owing 
to the interference effect from these two contributions. 

Among the seven modes (<j), tu, p° , ir°, rf\ fo)Ks, the first three are the states where final-state 
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TABLE V: Direct CP asymmetry parameter Af and the mixing-induced CP parameter A.Sj 
for various modes. The first and second theoretical errors correspond to the SD and LD ones, 
respectively (see the text for details). The foK$ channel is not included as we cannot make reliable 
estimate of FSI effects on this decay. 



Final State 




AS f 






.4,(96) 




SD 


SD+LD 


Expt 


SD 


SD±LD 


Expt 


<f>K s 

loK s 

P°K S 

V'K S 

j]K s 

TT°K S 


n 02+ 000 

u - U/: -0.04 

12+ 005 

-n nq +003 

u - uy -0.07 

n m +00 ° 

u - ui -0.04 

n n7 +0 - 02 

u - u ' -0.04 

06 +0 - 02 

u - UD -0.04 


n nq+o.oi+o.oi 

u - UD -0.04-0.01 
n m +0.02+0.02 
u - ul -0.04-0.01 
n n/1 +0.09+0.08 
U - U4 -0.10-0.11 

n nn +o.oo+o.oo 

u - uu -0.04-0.00 
n n7 +0.02+0.00 
u - u ' -0.05-0.00 
n ry+0.02+0.01 
u - u ^-0.03-0.01 


-0.38 ±0.20 

_n 17 +0.30 
Ul1 ' -0.32 

-0.30 ±0.11 

-0 3Q+ - 27 
u.oy_ 29 


1 4+0-3 

7 q+3.5 
' - J -2.6 

q n +2,2 
1 s+o- 4 

-fi 1+ 5 - 1 

-3.4+ 2 ;i 


9 fi +0.8+0.0 
0-0.4 
-i o 9+3.9+1.4 
io -^-2.8-1.4 

4fi fi+ 12 -9+3-9 
40 -°-13. 7-2.6 

9 -1+0.5+0.1 

z - x -0. 2-0.1 

q 7 +4.4+1.4 

°-< -1.8-2.4 

q 7+3.1 + 1.0 

°-' -1.7-0.4 


4± 17 

48 ±25 

4±8 
-8 ±14 



interactions may have a potentially large effect on the mixing-induced CP asymmetry Sf. In 
order to maximize the effects of FSIs on Sf, one should consider rescattering from charmless 
intermediate states that receive sizable tree contributions. Most of the intermediate states such 
as ~K°r]', K*~p + , ■ ■ •, in B decays are penguin dominated and hence will not affect Sf. For the 
decay B° — > 4>Ks, we have rescattering from K~ p + and K*~ir + . Because of the absence of 
the penguin chiral enhancement in B° — > K*~ir + and the large cancellation between 04 and ag 
penguin terms in B° — > K~ p + , it follows that the color-allowed tree contributions in these two 
modes are either comparable to or slightly smaller than the penguin effects. As for the luK$ mode, 
there is only one rescattering diagram, namely, B° — > K*~ir + — > K u, arising from the charmless 
intermediate states. (The rescattering diagram from K*~ p + is suppressed as elucidated before for 
the 4>Ks mode.) As a result, one will expect that final-state rescattering effect on Sf will be most 
prominent in B — > uKs, p K$ and <f>Kg- Indeed, we see from Table ITTT1 that FSI lowers SuK s by 
15% and enhances S p o Ks by 17% and S^Ks slightly. The theoretical predictions and experimental 
measurements for the differences between Sj D+LD and Sj/^k s , AS^ d+ld , are summarized in 
Table IV! It is evident that final-state interactions cannot induce large A5j in any of these modes. 

It is interesting to study the correlation between Af and Sf for the penguin dominated modes 
in the presence of FSIs. It follows from Eq. (jl.8|) that 

= cos 2(3 cot 5 f ps 0.95 cot S f , (3.21) 

for rf = (\ u Aj) / (\ c A C f) <C 1. This ratio is independent of |r,| and hence it is less sensitive to 
hadronic uncertainties. Therefore, it may provide a better test of the SM even in the presence of 
FSIs. Writing A u f ' c = \Ag'jj\e iS sD + iA u £^, we have 



\ X u^\A- SD co S 5^ + (\A u SD \ S ^ SD + Al D f 

\ r f\ = TT7T 1 

I* I v /|^ D cos^ D | 2 ±(|^ D |sin^ D + ^)2 

5 f = tan-^tan^ + ^ D ) - tan^tan^ + ^ D ), (3.22) 

\ A SD \ cos °SD I A SD I cos °SD 
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where the reality of Ald has been used, and Af and ASf can be obtained by using Eq. (|1.8|) . It 
is interesting to note that in the absence of LD contributions we have 1 5j D \ = \ Sg D — 5 C SD \ < ir/4 
for some typical SD (perturbative) strong phases and, consequently, we expect \ASf/Af \ > 1. 
This result generally does not hold in the presence of FSI. For example, in the case of l-A^pl <C 
\ A s'd\ ^ \ a ld\i [t is Possible to have \AS f /A f \<l. From Table |V| we obtain 

AS^kJA^Ks « -1-3 (1.5), ASuKs/AuKs « "0-08 (-1.7), AS p o Ks /A p o Ks « 0.08 (-1.1), 
iiSrfKslAfKs ~ -0-05 (0.6), ASiKg/AiKs ~ " 2 -0 (-1.1), ^S^kJA^o^ ~ 1-2 (-1-8), 

(3.23) 

where the SD contributions are shown in parentheses. For the foK$ mode, AS/ A ~ 3.0 but there 
is no reliable estimate of the FSI effect. From Eq. (|3.23f) we see that the relation \ASf/Af\ > 1 
indeed holds at short distances except for the rj'Ks mode where \$f\ ~ 65°. The sign flip of 
ASf/Af in the presence of LD rescattering is due to the sign switching of either Af (for 4>K$ and 
tt K s ) or AS f (for p°K s ). 



IV. CONCLUSIONS AND COMMENTS 



In the present work we have studied final-state rescattering effects on the decay rates and 
CP violation in the penguin-dominated decays B° — ► (4>,uj, p° ,rj' ,r],ir° , fo)Kg. Our main goal 
is to understand to what extent indications of possibly large deviations of the mixing-induced 
CP violation seen in above modes from sin 2(3 determined from B — > J/ipKs can be accounted for 
by final-state interactions. Our main results are as follows: 

1. We have applied the QCD factorization approach to study the short-distance contributions 
to the above-mentioned seven modes. There are consistently 2 to 3 a deviations between the 
central values of the QCDF predictions and the experimental data. 

2. The differences between the CP asymmetry Sf induced at short distances and the measured 
Sj/ipKs are summarized in Eq. (|2.18jl . The deviation of Sj D in the loKs and p°K$ modes 
from sin 2/3 is a 2 to 3 a effect owing to a large tree pollution. In contrast, tree pollution 
in rj'Ks is diluted by the QCD anomaly and/or final-state rescattering both of which carry 
negligible CP-odd phase. The long-distance effects on Sf are generally negligible except for 
the ujKg and p K$ modes where Sf is lowered by around 15% for the former and enhanced 
by the same percentage for the latter and ASfx^~pOA" s become consistent with zero within 
errors. 

3. Final state rescattering effects from charm intermediate states can account for the discrepancy 
between theory and experiment for the branching ratios of the modes uKs,r]'Ks,(f>Ks, and 
7r°Ks- Moreover, direct CP asymmetries in these modes are significantly affected; the signs 
of Af in the last two modes are flipped by final-state interactions. Direct CP asymmetries 
in the ujKs and p K$ channels are predicted to be A w k s ~ —0.13 and A p o Ks « 0.47, 
respectively, which should be tested experimentally. 
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4. For the fo(980)Ks mode, the short-distance contribution gives AS/A ~ 3.0, but at present 
we cannot make reliable estimates of FSI effects on this channel. 

5. Direct CP asymmetry in all the (b — ► s) penguin-dominated modes is rather small ( < a few 
%) except for ujKs and p°Ks- This strengthens the general expectation that experimental 
search for direct CP violation in b — > s modes may be a good way to look for possible effects 
of new CP-odd phases (see e.g. [3^). 

6. Since the mixing induced CP parameter Sf (actually ASf = —rjfSf — Sj/^x s ) an ^ the direct 
CP parameter Af are closely related, so are their theoretical uncertainties. Based on this 
study, it seems rather difficult to accommodate |A»Sy| > 0.10 within the SM, at least in the 
modes we study in this paper (except for foK$ which we cannot make reliable estimates). 

7. In particular, r/Kg and (to some degree) 4>Ks appear theoretically cleanest in our picture; 
i.e. for these modes the central value of ASf as well as the uncertainties on it are rather 



small. This also seems to be the case in QCDF [39j. Note also that the experimental errors 
on rj Kg are the smallest (see Table III) and its branching ratio is the largest, making it 
especially suitable for faster experimental progress in the near future |4jJ. 

8. The sign of ASf at short distances is found to be positive except for the channel p°Ks- 
After including final-state rescattering effects, the central values of ASf are positive for 
all the modes under consideration, but they tend to be rather small compared to the stated 
uncertainties so that it is difficult to make reliable statements on the sign at present. However, 
since the Sf and Af are strongly correlated, improved measurements could provide enough 
useful information that stronger statements on the sign could be made in the future. 
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